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Introduction 


Land  use  change  and  urban  development  are  known  to  effect  a  water¬ 
shed^  hydrologic  response  to  storm  rainfall.  Formerly  agricultural 
and  open  fields  become  partially  covered  with  an  impervious  surface  of 
rooftops,  driveways,  parking  lots,  and  roads,  connected  directly  to  the 
channel  by  storm  sewers.  This  transformation  of  land  cover  tends  to 
(1)  increase  the  total  volume  of  direct  runoff  by  decreasing  the 
opportunity  for  infiltration,  and  (2)  reduce  the  time  it  takes  runoff 
to  reach  the  channel  system  by  traveling  in  efficient  hydraulic  con¬ 
duits  rather  than  on  top  of  or  through  the  soil;  thus,  increasing  the 
peak  flow  in  the  stream.  As  one  result,  discharge  frequency  curves, 
derived  from  a  series  of  annual  peaks  that  occurred  under  varying  urban 
conditions,  will  be  inadequate  for  representing  the  present  probability 
of  flooding  in  the  basin. 

The  Chicago  District  and  North  Central  Division  were  concerned 
that  just  such  a  situation  exists  in  the  Dupage  River  Basin.  The 
District  has  contracted  with  an  A/E  firm  to  update  and  expand  a  HEC- 1 
model  of  the  basin  (Ref.  3).  The  A/E  was  required  to  calibrate  loss 
rates  against  available  discharge  frequency  curves.  Significant 
urbanization  has  taken  place  in  parts  of  the  basin  over  the  past  15 
years.  For  example,  between  1964  and  1975  rural  land  use  decreased  in 
one  subbasin  from  84  to  51  percent.  Annual  peaks  for  this  period 


represent  a  nonstationary  time-series  and  would  bias  both  discharge 
frequency  curves  and  HEC-1  model  parameters  calibrated  to  the  curves. 

The  Hydrologic  Engineering  Center  (HEC)  was  asked  in  July,  1977  by 
the  District  to  make  a  separate  study  of  the  basin  using  Hydrocomp's 
HSP,  a  continuous  hydrologic  simulation  computer  program.  HSP  can 
explicitly  account  for  spatial  and  temporal  variation  in  the  proportion 
of  impervious  area,  a  surrogate  measure  of  urban  development,  and  can 
generate  a  long  record  of  annual  peaks  from  observed  precipitation. 
Continuous  simulation  models  are  calibrated  only  to  observed  rainfall- 
runoff  events  and  not  probabilistic  estimates  of  annual  peak  dis¬ 
charge.  The  intended  products  from  HEC's  effort  were:  (1)  discharge 
frequency  curves  computed  from  simulated  peaks  and  representing  con¬ 
stant  present  conditions;  and  (2)  general  reconmendations ,  based  on  the 
modeling  experience,  concerning  the  value  and  requirements  of  continu¬ 
ous  hydrologic  simulation  in  urbanizing  basins. 

The  purpose  of  this  report  is  to  document  the  calibration  and 
application  of  an  HSP  model  for  the  DuPage  River  Basin,  and  to  evaluate 
HSP  as  a  tool  in  urban  hydrologic  analysis.  Topics  are  to  be  discussed 
in  the  following  order:  model  structure,  parameters,  and  data  require¬ 
ments;  a  description  of  the  study  basin,  parameter  calibration,  and 
land  use;  calibration  results,  long-record  simulation,  estimation  of 
discharge  frequency,  and  conclusions.  We  begin  with  a  brief  look  at 
how  the  model  works.  (A  complete  technical  description  of  HSP  can  be 
found  in  Reference  1.) 
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Model  Structure 


HSP  attempts  to  simulate  continuously  the  complete  hydrologic 
cycle  of  a  watershed.  Given  precipitation  and  potential  evapotrans- 
piration,  the  LANDS  module  accounts  for  the  following:  (1)  storage , 
as  interception,  upper  zone  (surface  or  near-surface),  lower  zone 
(sub-surface),  and  groundwater;  and  (2)  flow,  impervious  runoff, 
interflow,  overland  flow,  infiltration,  percolation,  base  flow  and 
actual  evapotranspiration.  Inflows  to  the  channel  system  are  routed 
downstream  in  the  CHANNELS  module  using  kinematic  wave  routing.  Figure 
1  is  a  schematic  of  HSP 9 8  major  components. 

The  model  has  the  additional  capability  to  simulate  the  complex 
snowfall-snowmelt  process.  Accumulation  of  precipitation  in  a  snowpack 
and  its  eventual  release  as  snowmelt  is  modeled  using  the  equations 
published  in  the  Corps  of  Engineers  MSnow  Hydrology”.  Melt  is  the 
combined  result  of  direct  solar  radiation,  convection,  condensation, 
rain  on  snow,  and  groundmelt.  Precipitation  and  temperature  decide 
when  snow  is  falling  and,  together  with  snow  density,  determine  depth 
and  equivalent  water  content  of  the  pack. 

Model  Parameters 


Sixteen  LANDS  parameters,  listed  and  defined  in  Table  l,  control 
the  capacity  of  storages  and  the  functional  relationship  between  flow 
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rates  and  storages.  For  some,  numerical  values  can  be  assigned 
directly  from  maps  or  other  data;  e.g.,  overland  flow  length  and  slope, 
groundwater  and  interflow  recession  constants,  percent  impervious  area, 
and  segment -to-gage  rainfall  ratio.  But  others,  including  the  param¬ 
eters  most  important  for  accurate  calibration  of  runoff  volume  and  peak 
flow  timing  and  magnitude,  are  not  physically  based  and  require  trial- 
and-error  adjustment.  Hydrocomp  suggests  "typical”  value  ranges  for 
upper  and  lower  zone  storage  nominal,  infiltration,  and  interflow. 
However,  several  runs  of  the  model  with  different  parameter  sets  and/or 
previous  HSP  experience  in  the  same  geographical  region  are  still 
necessary  for  accurate  calibration. 

Twelve  additional  parameters  are  involved  in  the  snowmelt  routine 
(Table  1).  Mean  elevation,  elevation  difference,  and  forest  cover  can 
be  measured  from  maps  but  the  other  nine  variables  are  difficult  to  set 
without  previous  experience  in  modeling  snowmelt  for  similar  meteoro¬ 
logical  conditions. 

The  channel  system  is  divided  into  reach  lengths  having  approxi¬ 
mate  constant  cross-sectional  geometry  and  hydraulic  roughness. 
Time-of-trave 1  at  bankfull  capacity  and  the  simulation's  computational 
time  step  are  also  considered  when  setting  up  the  CHANNELS  model.  For 
trapezoidal  channels,  each  reach  is  defined  by  9  parameters  (Table  2): 
channel  length  and  slope,  contributing  drainage  area,  top  and  bottom 
width,  depth,  flood  plain  slope,  and  Manning's  n  for  channel  and 
flood  plain. 
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Data  Requirements 

Hourly  precipitation  and  potential  evapotranspiration,  derived 
from  daily  pan  evaporation  data,  are  the  basic  input  data  required  by 
LANDS.  Continuous  daily  streamflow  and  hourly  hydrographs  for  selected 
storms  are  used  for  comparison  with  simulated  flows.  Meteorologic  data 
prerequisites  for  snowmelt  include  daily:  maximum  and  minimum  air 
temperature,  wind  movement,  dewpoint  temperature,  cloud  cover,  and 
incident  solar  radiation.  Point  source  discharges,  such  as  from 
municipal  or  industrial  sewage  treatment  plants,  can  contribute  a 
surprising  large  share  to  the  total  volume  of  annual  runoff.  To 
achieve  an  accurate  account  of  water  balance,  records  for  the  large 
point  sources  should  be  included  in  the  data  base. 

The  time-series  data  mentioned  above  provide  both  initial  input  to 
the  model  and  a  standard  to  evaluate  the  simulated  output. 

Supplemental  information  on  more  static  characteristics  of  the 
watershed  are  needed  for  estimating  parameter  values.  Land  use/land 
cover  help  determine:  impervious ness ,  forest  cover,  and  hydraulic 
roughness  for  overland  flow.  Slope,  length,  soil  permeability,  and 
channel  cross-section  geometry  define  several  parameters  in  LANDS  and 
CHANNELS.  When  basin  characteristics  change  with  time,  as  land  use 
does  in  urbanizing  basins,  multiple  measurements  may  be  necessary. 
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Description  of  Study  Basin 


r 


! 


The  DuPage  River  is  located  20  miles  west  of  Chicago  (Figure  2) 
and  has  a  drainage  area  of  324  square  miles  (measured  above  the  USGS 
stream  gaging  station  near  Shorewood,  Illinois).  The  East  and  West 
Branches  of  the  DuPage  River  join  below  the  town  of  Naperville  and 
together  account  for  nearly  two-thirds  of  the  basin  drainage.  The 
largest  concentrations  of  urban  land  are  along  the  East  Branch  and  in 
the  Cook  County  portion  of  the  West  Branch.  Below  their  confluence, 
the  mainstem  of  the  DuPage  River  collects  runoff  from  predominantly 
agricultural  and  rural  lands.  Continued  urban  development  has  been 
projected  for  both  East  and  West  Branches. 

It  was  originally  intended  for  HEC  to  study  the  entire  watershed 
and  to  determine  discharge  frequency  curves  at  all  of  the  basin's  six 
stream  gaging  stations.  However,  after  learning  how  to  use  HSP  and 
applying  the  model  to  one  DuPage  River  subbasin  it  was  decided  that 
time,  data,  and  funding  limitations  would  prevent  development  of  a 
complete  basin  model.  Rather,  one  subbasin  was  selected  for  full 
analysis:  West  Branch  DuPage  River  above  the  USGS  stream  gage  near 
West  Chicago.  The  selected  subbasin  (abbreviated  in  this  report  as 
"WBWC")  reflects,  in  miniature,  the  development  situation  facing  the 
entire  basin.  Recent  residential  construction  in  the  headwater 
villages  of  Hanover  Park  and  Schaumburg  has  changed  the  subbasin's 
hydrologic  behavior  and  could  signi f icant ly  effect  channel  flows 
through  lower  rural  areas.  Restricting  the  application  of  HSP  to  WBWC 
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should  not,  therefore,  interfere  with  the  study's  primary  objective:  to 
evaluate  a  continuous  hydrologic  simulation  model  under  changing  urban 
conditions . 


Calibration  of  Model  Parameters 


This  study  was  fortunate  in  being  able  to  acquire  a  data  base  for 
the  whole  basin  from  the  Northeastern  Illinois  Planning  Commission 
(NIPC).  A  product  of  a  previous  NIPC  water  quality  study  of  the  DuPage 
River,  the  data  base  contained  10  years  of  meteorologic ,  streamflow, 
and  point  source  data.  Table  3  lists  the  data  series  used  in  the 
calibration  of  the  WBWC  model.  Land  use  classification  at  three  points 
in  time,  maps,  reports,  and  channel  cross-sections  were  obtained  from 
NIPC  and  the  Chicago  District  and  were  very  useful  in  setting  up  the 
model  structure  and  estimating  model  parameters.  NIPC's  water  quality 
study  utilized  HSP's  LANDS,  CHANNELS,  AND  QUALITY  modules.  Initial 
values  of  many  variables  in  the  present  study  were  derived  from  NIPC's 
Hydrologic  Calibration  report  (Reference  2),  a  documentation  of  their 
own  application  of  LANDS  and  CHANNELS  to  the  DuPage  River. 

The  ten-year  data  base  was  divided  into  two  intervals,  1965-69  and 
1970-74,  and  analyzed  separately.  Model  parameters  were  calibrated 
using  the  earlier  period  data  and  a  four-step  iterative  procedure: 
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1.  Assign  new  parameters  values. 

2.  Run  simulation  on  5-year  calibration  data. 

3.  Compare  simulated  with  observed  streamflow. 

4.  Determine  which  LANDS  parameters  should  be  changed. 

Return  to  step  1. 

When  predicted  and  observed  runoff  were  sufficiently  close  the  model 
was  run  unchanged  on  the  remaining  5-year  verification  data.  If  this 
test  showed  major  systematic  bias  the  calibration  process  was 
restarted.  Otherwise,  the  model  was  accepted  as  a  valid  representation 
of  basin  hydrology. 

Initially,  four  land  surface  runoff  (LSRO)  segments  were  modeled: 
agricultural,  grassland,  forest,  and  impervious.  There  was  a  problem, 
however,  in  calibrating  with  such  a  detailed  land  use  breakdown.  At 
the  gaging  stations  only  total  runoff  from  all  surfaces  is  known  and 
not  the  individual  component  runoff  from  agricultural  land,  grassland, 
etc.  Differences  between  parameter  values  assigned  to  the  several  LSRO 
segments  cannot,  therefore,  be  justified  by  observed  data.  When 
parameters  were  set  to  what  appeared  to  be  reasonable  values,  it  was 
noticed  that  the  most  important  parameters  remained  the  same  for  all 
segments . 

Because  of  the  above  reasons,  a  simpler  two  segment  model  was 
adopted.  A  distinction  was  made  between  runoff  from  impervious 
surfaces  and  runoff  from  everthing  else,  referred  to  as  URBAN  and 
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RURAL,  respectively.  Since  only  impervious  runoff  is  permitted  in 
URBAN,  values  assigned  to  the  LANDS  parameters  for  this  segment  are 
arbitrary*.  In  contrast,  RURAL  allows  for  all  possible  runoff  except 
from  impervious  surfaces.  Its  final  LANDS  parameter  values,  including 
snowmelt,  are  given  in  Table  4. 

The  structure  of  CHANNELS  was  held  constant  for  all  computer 
runs.  The  WBWC  subbasin  was  split  into  three  channel  reaches  of  nearly 
equal  contributing  drainage  area.  Runoff  from  URBAN  and  RURAL  segments 
was  added  in  proportion  to  the  observed  land  use  of  each  reach. 
Cross-section  dimensions  and  the  other  CHANNELS  parameters  are  listed 
in  Table  4. 

Land  Use  Imperviousness 


Inventories  of  land  use/land  cover  reported  by  NIPC  for  1964, 

1970,  and  1975  provided  a  means  for  quantifying  change  in  the  amount  of 
impervious  surface  between  1965  and  1974.  First,  dissimilar  classifi¬ 
cations  were  collapsed  into  three  categories:  HI,  LO,  and  RURAL,  as 
shown  in  Table  5.  Then  land  use  was  translated  into  effective  imper- 
viou8ness  by  the  formula  IMP^CO.OLO  +  (0.2)HI.  An  attempt  was  made  to 
relate  imperviousness  and  land  use  by  actual  measurement  of  roof-tops, 
driveways,  streets,  etc.,  from  aerial  photographs.  This  approach, 
however,  produced  unrealistic  high  imperviousness  values.  The  above 

♦Except  for  A* 1 . 0  (  i .e . ,  100X  imperviousness)  and  the  snowmelt 
parameters . 
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simple  linear  equation  provided  reasonable  estimates  of  present 
imperviousness  and  would ,  at  a  minimum,  be  able  to  show  an  increase  in 
imperviousness  with  increasing  urban  development. 

Applying  the  imperviousness  relationship  to  1970  and  1975  land  use 
(Table  6.1),  percent  impervious  surface  (Table  6.2)  was  calculated  for 
the  total  subbasin  and  the  contributing  drainage  area  of  each  channel 
reach.  The  proportional  combination  of  land  surface  runoff  from  URBAN 
and  RURAL  segments  is  identical  to  the  impervious/pervious  breakdown  of 
Table  6.2. 

Calibration/Verification  Results 

The  HSP  model  was  evaluted  by  how  well  it  reproduced  runoff 
volumes  (annual,  monthly,  and  daily)  and  runoff  peak  discharges.  Even 
in  a  study,  such  as  the  present,  where  peak  flows  are  the  main  object 
of  interest,  it  is  necessary  to  be  able  to  generate  volumes  accurately 
over  long  periods  of  time.  This  is  because  of  the  significant  effect 
antecedent  storage  volumes  above  and  below  the  soil  surface  have  on  the 
volume  and  timing  of  storm  runoff.  Annual  and  monthly  volumes  are 
shown  in  Tables  7  and  8  for  the  calibration  and  verification  periods, 
respectively.  It  is  apparent  from  the  tables  that  March  to  April  were 
generally  low  in  simulated  volume  and  July  to  September  generally 
high.  Spring  runoff,  either  rain  plus  melting  snow  or  rain  on  frozen 
ground,  was  very  difficult  to  simulate  due  to  the  complex  nature  and 


data  uncertainties  of  the  snowmelt  process.  Runoff  from  summer 
thunderstorms  is  likewise  a  problem  because  of  the  hit-or-miss 
relationship  between  small  local  storm  cells  and  a  single  point 
observation  of  rainfall. 

Figures  3  and  4  are  accumulative,  or  mass  balance,  comparisons  of 
these  same  simulated  and  recorded  runoff  volumes.  They  show  that  after 
an  early  positive  error  corresponding  to  the  summer  storms  of  1963,  the 
mass  balance  remained  roughly  para  lied  to  the  “balance”  line  until  June 
1970.  Then  a  21  month  period  of  consecutive  high  volumes  began,  fol¬ 
lowed  by  2-1/2  years  of  lows.  Relative  to  the  calibration  period  the 
model  tested  poorly  for  volume  in  the  verification  years.  Partial 
explanation  of  this  change  in  performance  can  be  found  in  Figure  5.1 
where  annual  runoff  is  plotted  against  annual  precipitation.  Notice 
the  last  3  years  of  the  1970-74  period  had  significantly  higher 
runoff-to-precipitation  ratios  (indicated  by  the  slope  of  line  segments 
from  the  origin)  than  in  the  1965-69  period.  Using  published  data  for 
1975-77,  the  trend  for  an  increased  runoff-rainfall  ratio  is,  in 
general,  continued  to  be  observed.  It  is  reasonable  to  assume  the 
increase  in  runoff  per  unit  rainfall  is  caused  in  part  by  increased 
urbanization.  Figure  5.2,  a  plot  of  both  accumulative  annual 
precipitation  and  accumulative  annual  runoff  against  time,  supports 
this  hypothesis  by  contrasting  stationary  basin  rainfall  with  a  visably 
changed  incremental  runoff  volume.  The  break  in  slope  of  the  runoff 
curve  corresponds  directly  with  the  occurrence  of  large-scale 
residential  construction  in  the  upper  region  of  the  subbasin. 
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As  mentioned  before  in  the  land  use  discussion,  a  change  was  made 
in  the  proportion  of  impervious  surface  from  2.62  in  1965-69  to  5.5%  in 
1970-74.  It  may  be  that  a  larger  increase  was  necessary  to  generate 
higher  runoff  volume  for  the  latter  period.  To  be  sure  a  detailed 
examination  of  several  individual  storm  volumes  should  be  made.  Data 
limitations  prevented  such  an  analysis  from  being  included  in  the 
present  study. 

Mean  daily  flows  can  be  considered  as  both  average  flow  rate  (cfs) 
and  daily  runoff  volume  (cfsd).  Continuous  hydrographs  for  water  years 
1965-74,  Figures  6.1  to  6.10,  present  a  graphical  comparison  of  simu¬ 
lated  and  observed  daily  flows.  Plots  such  as  these  were  routinely 
generated  during  calibration-  They  provided  an  opportunity  to  see  not 
only  peaks  and  total  volumes  but  also  the  relative  contribution  of 
component  flows;  e.g.,  surface  runoff,  interflow,  and  groundwater.  The 
average  of  the  simple  correlation  coefficients  (r),  calculated  sepa¬ 
rately  for  each  water  year  on  Figures  6.1  -  6.10,  was  0.82. 

A  primary  study  objective  was  the  generation  of  long  record 
synthetic  peak  discharge.  Consequently,  a  great  deal  of  attention  was 
given  to  the  accurate  reproduction  of  large  peak  flows.  Tables  9  and 
10  list  observed  and  simulated  peak  discharge  for  events  corresponding 
to  the  USGS  published  partial-duration  series*.  Mean  error  for  the 

♦Not  included  in  the  tables  are  events  for  which  either  (1)  the  simu¬ 
lated  peak  was  not  printed  on  the  computer  output,  or  (2)  the  observed 
discharge  was  below  the  specified  base  of  the  partial-duration  series. 
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calibration  period  was  much  smaller  than  for  the  verification  period 
(+2. OX  vs  -22. 6Z),  with  both  periods  having  about  equal  "spread"  in  the 
distribution  of  error.  Simulated  peaks  for  water  years  1971-74  were 
consistently  low.  An  explanation  identical  to  that  provided  for  the 
occurrence  of  low  flow  volumes  applies  to  the  present  case;  i.e.,  a 
combined  result  of  unusually  high  runoff  (for  a  given  precipitation) 
and  inadequate  increase  in  impervious  surface. 

Long-Record  Simulation 

When  the  incremental  improvement  in  model  fit  became  small  for 
successive  changes  in  parameter  values,  the  calibration  process  was 
stopped.  The  final  set  of  parameter  values  and  the  percent  imperv¬ 
iousness  used  in  the  verification  runs  (i.e.,  "present  conditions") 
together  make  the  HSP  model  a  mechanism  for  translating  observed 
precipitation  into  streamflow  at  the  subbasin  outlet.  Having  a 
calibrated  model,  the  next  task  was  to  assemble  a  long  continuous 
record  of  the  necessary  meteorologic  variables. 

Both  the  length  and  accuracy  of  a  simulation  are  frequently 
limited  by  the  data  available,  or  rather,  unavailable.  There  were  very 
few  hourly  recording  stations  in  existence  25  years  ago.  And  even  if 
collected,  precipitation  data  did  not  begin  to  be  systematically  coded 
in  computer-compat ible  format  by  the  National  Weather  Service  until 
1948.  Depending  upon  local  spatial  variation  in  precipitation,  the 
quality  of  a  simulation  may  be  adversely  affected  by  the  distance  to 
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and  location  of  the  rain  gage(s).  These  same  factors  apply  to  the 
other  required  meteorologic  data. 

The  longest  hourly  precipitation  record  in  the  vicinity  of  the 
basin  is  station  1582,  Chicago  WB  City,  located  in  downtown  Chicago 
near  the  lakefront.  Data  has  been  collected  here  since  July  1899,  but 
available  on  magnetic  tape  only  from  September  1948.  Therefore,  the 
period  of  long-record  simulation  was  determined  by  default  to  be  water 
years  1949-64*.  Pan  evaporation  data  for  this  same  period  was  not 
readily  available.  Instead  potential  evapotranspirat ion  at  O’Hare 
Airport  for  1965-74,  as  used  in  the  calibration/verification  phase,  was 
duplicated  twice  to  create  an  artificial  potential  evapotranspirat ion 
record  for  1949-64.  A  similar  procedure  was  followed  to  construct 
radiation  and  wind  speed  data  files  for  the  long-record  years.  Minimum 
and  maximum  daily  temperatures  were  measured  at  the  Chicago  WB  City 
station  and  were  used  in  the  snowmelt  calculations.  Table  11 
summarizes  the  long  record-simulation  input  data. 


Table  12  contains  the  simulated  monthly  and  annual  flow  volumes 
(ail  generated  for  constant  ’’present  conditions'’)  for  the  following 
time  intervals:  (l)  long-record  period,  1949-64; 


♦Because  the  Chicago  WB  City  station  was  discontinued  in  November 
1964,  a  double  mass  plot  of  accumulative  precipitation  for  Chicago  WB 
City  vs.  Roselle  could  not  be  made.  Data  was  available  for  two  years, 
1967-68,  for  the  replacement  station  Chicago  WB  City  #2.  However,  due 
to  the  short  concurrent  period,  a  double  mass  plot  against  the  Roselle 
station  was  not  conclusive.  Without  additional  informat ion,  we  were 
forced  to  assume  that  the  long-record  Chicago  WB  City  data  is  repre¬ 
sentative  of  the  true  basin  precipi tat ion  for  the  period  1949-64. 


14 


(2)  calibration  period,  1965-69,  rerun  under  present  land  use  condi¬ 
tions;  and  (3)  verification  period,  1970-74,  identical  to  values 
previously  reported  in  Table  8.  Table  13  lists  both  simulated  annual 
peak  discharge  for  these  same  periods  and  the  corresponding  published 
values  of  observed  peak  flow. 

Discharge  Frequency 

Discharge  frequency  curves  were  computed  from  the  simulated  and 
observed  annual  peak  series  of  Table  13,  and  plotted  as  Figure  7.1. 

The  curves  represent  the  probability  distribution  of  annual  peak 
discharge  for  two  distinct  populations:  the  real  (observed)  world  and 
the  model  (simulated)  world.  Although  the  two  curves  appear  to  be 
different  (which  is  what  we  would  expect  them  to  be,  as  discussed 
later),  it  will  be  useful  to  measure  how  significant  their  difference 
really  is. 

Discharge  frequency  curves  are  defined  by  assuming  a  distribution 
(in  our  case,  log-Pearson  Type  III)  and  estimating  the  necessary 
parameters  (mean,  standard  deviation,  and  skew).  A  generalized  skew 
value  of  0.0  was  adopted  for  the  present  study.  Hypothesis  tests  were 
performed  on  the  remaining  parameters  to  determine,  in  probabilistic 
terras,  if  the  differences  in  sample  statistics  were  significant*. 

Using  a  t-test  to  examine  the  hypothesis  H:  jj  ^  sy  ^  (that  is, 

♦Details  of  the  tests  are  given  in  Appendix  A. 
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equivalence  of  means),  it  was  found  that  at  a  significance  level 

a  3  0.40  the  hypothesis  could  not  be  rejected.  This  implies  that  even 

with  the  test  structured  so  that  there  was  a  40%  chance  of  rejecting  a 

true  hypothesis  (Type  I  error)  the  equivalence  of  means  was  not 

2  2 

rejected.  An  F-test  of  the  hypothesis  H:  a  ^  =  a  2  (e9ual 

variance)  showed  that  at  a  significance  level  a  =  0.50  the 
null  hypothesis  once  again  could  not  be  rejected.  From  a  statistical 
viewpoint,  therefore,  the  difference  between  the  two  sample  means  and 
the  difference  between  the  two  sample  variances  are  not  significant. 

Despite  the  statistical  comparison,  there  are  several  reasons  why 
we  should  expect  the  curves  to  be  different. 

(1)  The  observed  frequency  curve  is  based  on  a  series  of  annual 
peaks  that  occurred  under  varying  degrees  of  urbanization 
whereas  the  simulation  series,  and  the  frequency  curve 
associated  with  it,  was  deliberately  intended  to  represent 
constant  Mpresentn  urban  conditions. 

(2)  The  simulation  record  consisted  of  twenty-six  annual  peaks 
for  the  period  1949-74.  This  is  a  longer  and  different  set 
of  years  than  the  observed  record's  eighteen  peaks  from  the 
period  1961-78. 

(3)  The  model  is  a  simplification  of  the  real  world's  hydrologic 
and  hydraulic  processes.  It  uses  discrete  point  measurements 
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of  the  basic  input,  precipitation,  and  attempts  to  generalize 
this  information  both  in  time  and  space.  It  must,  therefore, 
generate  peak  flows  that  do  not  agree  exactly  with  the  actual 
storm  events. 

(4)  As  described  previously  in  the  section  on  long-record 
simulation,  both  a  different  recording  station  (e.g., 
precipitation  and  temperature)  and  an  artificial  data  record 
(e.g.,  potential  evapotranspiration  and  solar  radiation)  were 
used  to  generate  streamflow  for  the  period  1949-64.  The 
magnitude  of  the  error  introduced  into  the  simulation  by  the 
non-representativeness  of  these  data  is  unknown,  but  such 
error  surely  exists. 

An  attempt  was  made  to  examine  separately  some  of  the  possible 
explanations  of  why  the  curves  of  Figure  7.1  should  be  different.  In 
Figure  7.2  is  a  plot  of  the  simulated  and  observed  frequency  curves  for 
the  same  ten  years,  1965-74.  Recall  that  this  was  the  calibration/ 
verification  period  for  the  simulation  model,  and  that  during  this 
period  the  model's  percent  imperviousness  was  changed  from  2.6%  in 
1965-69  to  5.5%  in  1970-74.  This  transition  was  made  to  try  to  imitate 
the  actual  change  in  urbanization  that  was  occurring  in  the  watershed. 
Hypothetically,  the  two  curves  of  Figure  7.2  should  represent  the  same 
changing  urban  conditions  and  differ  only  in  the  degree  to  which  the 
model  was  unable  to  imitate  the  real  world. 


As  before,  hypothesis  tests  were  made  on  the  equivalence  of  means 

and  the  equivalence  of  variances.  It  was  found  that  the  hypothesis 

H:  y  ^  *  y  £  could  not  be  rejected  at  the  a  =  0.50  significance 

2  2 

level,  and  the  hypothesis  could  not  be  rejected  at  the 

a  =0.20  significance  level.  The  conclusion  that  should  be  made  is 
the  following:  the  small  sample  size  (n  =  10)  does  not  permit  the 
inference  (from  the  observed  differences  in  sample  means  and  variances) 
that  the  corresponding  population  parameters  are  significantly 
di ff erent . 

However,  visual  examination  of  Figure  7.2  shows  that  at  the  tails 
of  the  distribution  (low  and  high  exceedance  probabilities)  the  curves 
do  depart.  For  example,  at  the  IX  event  the  simulated  frequency  curve 
is  230  cfs  higher  than  the  observed  frequency  curve  (1310  vs.  1080). 

To  put  this  magnitude  in  perspective,  230  cfs  i9  more  than  twice  the 
difference  between  the  1%  and  2%  flood  events  of  the  observed  curve 
(1080  and  982,  respectively).  The  plotted  annual  peaks  for  this  time 
period  reveals  a  general  trend:  the  three  largest  simulated  peaks 
exceed  or  are  nearly  equal  to  the  observed  peaks,  whereas  the  remaining 
seven  simulated  peaks  are  all  lower  than  the  comparable  observed 
peaks.  Rather  than  random  variation  about  the  observed  data  there 
appears  to  be  a  systematic  bias  in  the  calibrated  model  of 
over s imulat ing  large  events  while  undersimulating  small  ones. 

Figure  7.3  illustrates  the  effect  of  the  longer  simulation 
record.  Discharge  frequency  curves  were  computed  from  the  simulated 
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peaks  (at  constant  present  urban  conditions)  for  two  overlapping 

periods:  1949-74  (26  years)  and  1961-74  (14  years).  Both  hypotheses, 

2  2 

H:  P  ^  *  Vij  an<*  H:  °i  ®  °2  f  were  not  rejected  at 

the  a  ■  0.50  significance  level.  This  result  should  have  been  expected 
as  the  two  curves  were  derived  from  non-d is joint  samples  of  the  same 
population.  Notice  in  Figure  7.3  the  difference  between  the  simulated 
frequency  curves  based  on  different  lengths  of  record  is  small  relative 
to  the  difference  between  simulated  and  observed  curves  for  the  same 
period,  1961-74. 

The  simulated  frequency  curve  of  Figure  7.1  was  intended  to  show 
(by  means  of  the  simulated  peaks)  what  would  have  been  the  actual 
frequency  curve  if  urbanization  were  held  constant  at  its  present 
(1974)  level  and  annual  peak  data  were  available  for  the  longer  1949-74 
period.  That  hypothetical  present-day  (1974)  frequency  curve  is 
contrasted  in  Figure  7.1  with  the  curve  derived  from  historically 
observed  peak  streamflows.  From  the  above  analysis  it  seems  unlikely 
that  changing  urban  conditions  and/or  a  longer  record  can  explain  all 
the  differences  between  the  two  curves.  Differences  due  to  these 
causes  have  very  likely  been  overshadowed  by  the  model,  as  calibrated, 
generating  annual  peak  events  that  contain  the  systematic  bias 
previously  discussed. 
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Conclusions 


The  application  of  HSP  to  urban  watersheds  for  the  purpose  of 
determining  discharge  frequency  can  be  evaluated  by  the  following 
criteria; 

(1)  How  accurate  was  the  model's  simulation?  Accuracy  can  be 
judged  at  four  levels:  monthly  and  annual  volumes,  daily 
flows,  peak  discharge,  and  discharge  frequency. 

(a)  The  average  absolute  monthly  volume  error  (Tables  7  and 
8)  was  32.12  and  28.12  for  the  calibration  and 
verification  periods,  respectively.  A  seasonal  pattern 
of  underestimating  volumes  in  the  spring  and 
overestimating  volumes  in  the  late  summer  tended  to 
cancel  each  other,  as  evidenced  by  the  smaller  absolute 
yearly  volume  errors  of  11.12  and  16.12  for  calibration 
and  verification,  respectively. 

(b)  Comparing  daily  flows  will  show  how  the  model  performed 
in  distributing  runoff  volume  into  daily  increments.  The 
correlation  of  mean  daily  flows  within  each  of  the  ten 
calibration/verification  water  years,  Figures  6.1  to 
6.10,  ranged  from  0.75  to  0.92,  and  averaged  0.82. 
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(c)  Peak  discharges  for  27  events  in  the  calibration  period 
and  35  events  in  the  verification  period  were  available 
for  comparison*  The  5-year  calibration  group  (Table  9) 
had  an  average  absolute  error  of  26.0%,  High  and  low 
errors  tended  to  compensate  for  each  other,  producing  a 
simple  average  error  of  2.0%.  Peaks  from  the  5-year 
verification  period  (Table  10)  had  an  average  absolute 
error  of  36.0%  and  a  simple  average  error  of  -22.6%, 
demonstrating  a  consistent  low  bias  for  nearly  all  storm 
events  in  this  period.  Possible  reasons  for  this 
behavior  of  the  model  were  previously  discussed. 

(d)  A  series  of  statistical  tests  on  the  discharge  frequency 
curve  parameters  (Figure  7.1)  concluded  that  differences 
between  sample  statistics  were  not  significant. 
Examination  of  the  potential  reasons  as  to  why  the  curves 
should  be  (and  visually  appear  to  be)  different 
demonstrated  that  error  in  the  calibrated  model  was  the 
primary  cause  of  this  difference. 

(2)  Can  the  model  account  for  changing  urbanization? 

Theoretically  yes  -  the  model  can  be  structured  (e.g., 
separate  pervious  and  impervious  land  surface  runoff 
segments)  to  conveniently  account  for  changing  urban 
conditions.  The  problem  arises,  as  is  the  case  in  all  urban 
hydrologic  models,  in  going  from  known  land  use/land  cover  to 
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the  amount  of  effective  impervious  surface.  Other  changes 
associated  with  urban  areas,  e.g.,  channel  improvements  or 
flood  control  structures,  can  be  included  in  the  model. 

(3)  Is  the  model  easy  to  use? 

To  operate  -  yes,  to  calibrate  -  no.  Model  structure, 
parameter  values,  and  input  data  to  be  used  and  output  data 
to  be  saved  are  all  specified  by  simple  program  statements. 
The  storage,  retrieval  and  manipulation  of  large  time-series 
data  structures,  and  the  selection  of  graphical  and 
statistical  operations  are  defined  with  a  set  of  program 
commands  written  in  natural  language  with  a  logical 
well-defined  syntax.  Without  previous  experience,  however, 
calibration  to  a  specific  study  basin  can  be  a  difficult  and 
t ime-consuming  task.  There  are  sixteen  parameters  in  LANDS 
to  manipulate  and  with  snowmelt  almost  twice  that  many. 
Several  LANDS  parameters  initially  assigned  on  the  basis  of 
maps  or  other  data  were  found  later  to  require  changing  to 
obtain  an  acceptable  fit.  The  most  important  variables 
(UZSN,  LZSN,  INTERFLOW,  INFILTRATION)  are  not  physically 
derived.  Without  having  experience  with  the  model  in  a 
similar  hydrologic  area  it  is  very  hard  to  know  what  values 
to  start  with  for  these  parameters,  or  how  much  to  change 
them  on  successive  computer  runs. 
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(4)  What  are  the  data  requirements? 

Extensive,  both  in  quantity  of  data  and  hours  of  labor 
necessary  for  preprocessing.  Hourly  precipitation  and  daily 
evaporation  data  are  required  for  input  and  daily  streamflow 
to  evaluate  simulation  output.  Hourly  precipitation  prior  to 
1948  is  seldom  available  in  computer  format  and  the  enormous 
job  of  locating,  coding,  cleaning,  and  keypunching  could  be 
restrictive  in  some  studies.  The  snowmelt  routine  requires 
at  least  semi-daily  temperature;  if  not  available,  daily 
radiation,  wind  speed,  and  dewpoint  can  be  estimated. 
Temperature  is  usually  measured  at  many  stations  but  the 
other  meteorologic  variables,  including  evaporation,  can  be 
difficult  to  obtain. 

(5)  Can  runoff  from  several  different  land  surfaces  be  modeled? 

In  theory  -  yes.  The  model  is  structured  so  that  runoff 
caused  by  any  possible  combination  of  hydrologic  and 
meteorologic  factors  can  be  generated  separately;  e.g.,  from 
different  rain  gages  for  forest,  pasture,  residential,  and 
agricultural  land  cover.  The  problem  in  using  this 
capability  is  not  knowing  what  the  individual  runoff  from 
each  segment  should  be,  only  their  aggregate  contribution  at 
the  streamgage.  Parameter  modifications  for  multi-land  use 
watersheds  must  therefore  be  based  on  rational  thought  rather 
than  empirical  evidence. 
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(6)  What  additional  benefits  are  there  in  continuous  simulation? 
One  benefit  is  a  greater  understanding  of  and  appreciation 
for  the  complexity  of  a  watershed’s  hydrologic  cycle. 

Relative  contributions  of  surface  runoff,  interflow,  and 
groundwater  are  noticed  as  are  evaporation,  infiltration,  and 
percolation.  The  seasonal  pattern  of  rainfall-runoff  and 
their  variations  in  volume  from  year-to-year  are  made 
explicit.  Perhaps  having  improved  knowledge  of  a  basin's 
hydrologic  behavior  will  permit  a  greater  number  and/or 
variety  of  water  resource  management  alternatives  to  be 
considered. 

(7)  How  long  does  it  take  and  what  does  it  cost  to  develop  an  HSP 
mode 1 ? 

Starting  with  a  nearly  complete  hydrologic/meteorologic  data 
base,  it  took  one  person  6  months  to  assemble  and  analyze 
additional  data,  and  to  learn  how  to  use  the  model.  Another 
6  months  were  spent  in  calibration  and  long-record  simula¬ 
tion.  Time  gained  by  experience  (e.g.,  shortening  of 
calibration  time)  would  probably  be  canceled  by  what  would 
have  been  required  to  construct  the  original  data  base.  A 
realistic  estimate  for  studying  a  basin  similar  to  the  DuPage 
River  would  be  9  to  12  person-months.  The  cost  for  a  single 
LANDS  simulation  of  the  5-year  calibration  period  was  about 
$30.  Total  computer  cost  should  probably  be  less  than 
$2,000.  Although  these  figures  apply  to  the  modeling  of  a 
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2 

28.5  mi  basin,  it  is  expected  that  cost  would  increase 
less  than  proportionately  when  modeling  a  larger  watershed. 

(8)  Are  continuous  simulation  models  well  suited  for  studies 
involving  only  the  estimation  of  discharge  frequency  in 
urbanizing  basins? 

They  do  the  job,  but  considering  the  amount  of  work  (and 
data)  required,  it  is  doubtful  that  they  are  the  best  tool 
for  such  limited  studies.  More  comprehensive  study 
objectives  (for  example,  water  supply,  flood  control,  or 
water  quality,  which  explicitly  consider  runoff  volume  and 
water  balance)  would  take  advantage  of  the  continuous 
simulation’s  detailed  analysis.  Discharge  frequency  under 
changing  urban  conditions  is  a  problem  that  could  be  handled 
by  simpler,  quicker,  less  costly  approaches  requiring  much 
less  data;  e.g.,  design  storms  or  several  historical  events 
used  as  input  to  a  single-event  model,  or  a  continuous  model 
with  a  less  complex  soil-moisture  accounting  algorithm. 

In  summary,  an  HSP  continuous  hydrologic  simulation  model  of  the 
West  Branch  DuPage  River  above  West  Chicago  was  calibrated  and 
verified,  long-record  precipitation  and  other  meteorologic  data  were 
assembled  and  provided  as  input  to  the  model,  and  resulting  annual  peak 
discharges  used  to  estimate  discharge  frequency  at  the  gage  site  for 
constant  present  conditions.  Considering  the  major  strengths  and 
weaknesses  identified  above,  the  following  conclusion  can  be  made. 
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Continuous  hydrologic  simulation  provides  a  viable  alternative  as  a 
method  of  analysis  for  urban  hydrologic  systems.  It  requires  large 
amounts  of  data  and  significant  labor  to  calibrate  an  HSP  model  to  a 
specific  basin*  For  studies  which  warrant  such  a  detailed  definition 
of  the  hydrologic  process,  the  model  provides  a  rational  representation 
of  basin  hydrology  that  could  be  applied  to  the  investigation  of 
complex  water  resource  problems. 
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LAND 

K1 

A 

EPXM 

UZSN 

LZSN 

K3 

K24L 

K24EL 

INFILTRATION 

INTERFLOW 

L 

SS 

NN 

IRC 

KV 

KK24 


SNOW 

RADCON 

CONDS-CONV 

SCF 

ELD  IF 

IDNS 

F 

DGM 

WC 

MPACK 

EVAPSNOW 

MELEV 

TSNOW 


*  Adopted  from 


TABLE  1* 

DUPAGE  RIVER  BASIN 
MODEL  PARAMETER  DEFINITIONS 
LANDS 


Ratio  of  average  segment  rainfall  to  average  gage 
Impervious  area  (fraction) 

Interception  storage  (maximum  value) 

Nominal  upper  zone  soil  moisture  storage 
Nominal  lower  zone  soil  moisture  storage 
Actual  evaporation  rate  parameter 
Seepage  to  'deep*  groundwater 
Evaporation  from  perched  groundwater 
Infiltration 
Interflow 

Length  of  overland  flow 
Overland  flow  slope  (ft/ft) 

Manning* s  N  for  overland  flow 
Daily  interflow  recession  rate 
Groundwater  recession,  variable  rate 
Groundwater  recession,  constant  rate 


Radiation  melt  parameter 
Convection  melt  parameter 
Snow  correction  factor  to  gage  record 
Elevation  difference  (gage  to  segment) 

Initial  density  of  new  snow 
Forest  cover 

Daily  ground  melt  (inches) 

Water  content  of  snowpack  maximum 
Snowpack  at  complete  areal  coverage 
Snow  evaporation  parameter 
Mean  watershed  segment  elevation  (ft) 

Upper  limit  of  temperature  at  which  precipitation  is 


reference  2. 


snow 
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REACH 

LIKE 

TYPE 


TO# 

SEG# 

LEN 

AREA 

UPSTR 

DNSTR 

W1 


W2 

DEPTH 

S-FP 

N-CH 

N-FP 


*  Adopted 


TABLE  2* 

DUPAGE  RIVER 

MODEL  PARAMETER  DEFINITIONS 
CHANNEL 


Reach  Number 

Reach  number  that  has  an  identical  cross  section 
The  type  of  channel: 

RECT:  Trapezoidal  channel  cross  section 

CIRC:  Circular  Conduit 
IMAG:  Feeder  reach  without  routing 

DAM  :  Reservoir 

Reach  number  to  which  the  reach  is  tributary 
Land  surface  segment  that  contains  the  reach 
Length  of  the  reach  in  miles 

Local  area  tributary  to  the  reach  in  sq.  miles 

Upstream  channel  bottom  elevation  in  the  reach 

Downstream  channel  bottom  elevation  in  the  reach 

Incised  channel  bottom  width  in  feet  for  trapezoidal 
channels,  or  the  diameter  in  inches  for  circular 
channels 

Incised  channel  top  width  in  feet  for  trapezoidal 
channels  or  Manning's  N  for  circular  channels 

Incised  channel  depth  in  feet 

Transverse  slope  of  the  flood  plain  in  feet  per  foot 
Manning's  N  for  the  incised  channel 
Manning’s  N  for  the  flood  plain 


from  reference  2. 
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Data 

Precipitation 

Evaporation 

Temperature 

Wind 

Cloud  Cover 
Dewpoint 
Radiation 
Streamflow 

Point  Source 
Point  Source 


TABLE  3 


Calibration  Data 

Interval 

Location 

Hourly 

Rose  lie 

Daily 

Chicago  O' Hare 

Semi -Daily 

Chicago  O'Hare 

Daily 

Chicago  O' Hare 

Daily 

Chicago  O'Hare 

Daily 

Chicago  O' Hare 

Daily 

Chicago  O'Hare 

Daily 

West  Branch  DuPage  River 
near  West  Chicago 

Daily 

Hanover  #1  STP 

Daily 

Hanover  Park  STP 

Point  Source 


Old  Bartlett  STP 


TABLE  4 


TABLE  5 


Land  Use/Land  Cover  Classifications 


1964 

NIPC 

1970 

NIPC 

1975 

LANDSAT 

PRESENT 

STUDY 

Conmerc  ial 
Industiral 
Mining 

Manufacturing 

Trade 

Services -Private 
Warehouse 

Shopping  Centers 
Hotel/Motel 

Parking 

Mining 

High-Intensity 
Developed  Land 

HI 

Residential 
Streets 
Transpo/Com/ 
Utilities 
Public  Bldgs 

Res-Single  Fam 
Res-Milti  Fam 

Mobile  Homes 

Streets 

Trans po /Com/Util 
Airports 

Railroads 

Services-Inst 

Public  Bldgs 

Military 

Low-Intensity 
Developed  Land 

LO 

1 

Public  Open 
Space 

Agr ic/ Vacnt 

Water 

Public  Open  Space 
Entertainment 
Assembly 

Cemetary 

Vacant-Under  Devlp 

Vacant/ Agr ic/ For st 

Water 

1 

Grass 

Trees 

Soybeans 

Corn 

Other  Crops 

Open  Space 

Water 

RURAL 

i 

1 

i 

1 


Table  6.1 

Land  Use  Distribution  (Percent) 


Reach 

# 

DA 

(mi2) 

1970 

i 

1975  j 

HI 

LO 

RURAL 

HI 

LO 

RURAL 

5 

10.2 

1.8 

39.1 

59.1 

9.0 

59.9 

31.1 

10 

9.0 

1.1 

15.6 

83.3 

7.0 

49.6 

43.4 

15 

9.3 

1.4 

10.7 

87.9 

1.2 

18.4 

80.4 

WBWC 

28.5 

1.5 

22.7 

75.8 

5.9 

43.5 

50.6 

Table  6.2 


Impervious  Surfac  Distribution  (Percent) 


Reach 

# 

DA 

(mi2) 

1970 

1975 

Impervious  1 

Pervious2 

Impervious^ 

Pervious2 

5 

10.2 

4.3 

95.7 

7.8 

92.2 

9.0 

1.8 

98.2 

6.4 

93.6 

15 

9.3 

.1,4 

98.6 

2.1 

97.9 

WBWC 

28.5 

2.6 

97.4 

5.5 

94.5 

1.  Impervious  =  (.1)  LO  +  (.2) HI 

2.  Pervious  *  100  -  (Impervious) 
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TABLE  9 


1 


PEAK  DISCHARGE  COMPARISON 
CALIBRATION  PERIOD  1965-69 


WATER 

Peak  Discharge  (cfs) 

DIFFERENCE 

DIFFERENCE 

YEAR 

DATE 

RECORDED 

SIMULATED 

(  C  f  8  ) 

(Z) 

65 

1/2 

151 

174 

23 

15.2 

1/23 

116 

164 

48 

41.4 

2/7 

250 

207 

-43 

-17.2 

4/1 

180 

147 

-33 

-18.3 

4/6 

196 

215 

19 

9.7 

4/15 

113 

104 

-9 

-8.0 

4/25 

105 

119 

14 

13.3 

66 

12/25 

151 

18 

13.5 

2/10 

152 

-219 

-59.0 

4/21 

114 

8 

7.5 

4/24 

124 

140 

16 

12.9 

4/28 

144 

315 

171 

118.8 

5/12 

537 

474 

-63 

-11.7 

i 

1 

67 

4/1 

461 

-206 

! 

-44.7 

4/30 

120 

43 

35.8 

6/10 

805 

912 

107 

13.3 

6/17 

201 

114 

-87 

-43.3 

6/28 

133 

112 

-21 

i 

-15.8 

1 

68 

8/17 

340 

651 

i 

311 

91.5 

1 

69 

11/29 

111 

115 

4 

3.6 

1/23 

223 

203 

-20 

-9.0 

1/30 

142 

138 

-4 

-2.8 

4/5 

183 

130 

-53 

;  -29.0 

4/18 

192 

111 

-81 

-42.2 

6/9 

399 

358 

-41 

-10.3 

7/1 

100 

102 

2 

;  2.0 

7/20 

i 

177 

i 

154 

-23 

-13.0 

Average  Error  2.0 


Average  Absolute  Error  26.0 
Standard  Error  38.2 
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TABLE  10 

PEAK  DISCHARGE  COMPARISON 
VERIFICATION  PERIOD  1970-74 


WATER 

Peak  Discharge  (cfs) 

DIFFERENCE 

DIFFERENCE 

YEAR 

DATE 

RECORDED 

SIMULATED 

(cfs) 

(Z) 

70 

mam 

237 

186 

-51 

-21.5 

.■in  !:M 

205 

286 

81 

i  39.5 

521 

438 

-83 

-15.9 

317 

241 

-76 

-24.0 

180 

305 

125 

69.4 

9/7 

186 

400 

214 

115.1 

9/25 

207 

L.  . 

199 

-80 

-3.9 

71 

2/20 

333 

120 

,  -213 

;  -64 . 0 

72 

12/15 

i 

344 

259 

-85 

;  -24.7 

3/13 

209 

118 

-91 

-43.5 

A/17 

404 

229 

-175 

-43.3 

4/22 

324 

135 

-189 

-58.3 

6/15 

562 

231 

-331 

-58.9 

6/20 

555 

292 

-263 

|  -47.4 

7/18 

380 

130 

|  -250 

I  -65.8 

i  8/26 

715 

687 

-28 

-3.9 

9/14 

438 

288 

-150 

-34.2 

9/18 

448 

377 

-71 

-15.8 

9/29 

378 

290 

-88 

-23.3 

73 

10/12 

195 

153 

-42 

-21.5 

10/23 

2  59 

138 

-121 

-46.7 

12/30 

535 

444 

-91 

-17.0 

4-22 

383 

300 

-83 

-21.7 

5/28 

154 

103 

-51 

-33.1 

6/17 

158 

i 

-  .  .  i 

177 

19 

12,0 

74  j 

10/13 

1 

234  | 

157 

-77 

-32.9 

12/5 

234 

125 

-109 

-46.6 

l 

1/21 

244 

227 

-17 

-7.0 

1/27 

522 

314 

-208 

-39.8 

j 

3/5  J 

220 

187 

-33 

-15.0 

1 

4/14 

167 

111 

-56  ! 

-33.5 

5/17 

472 

182 

-290 

-61.4 

5/22 

376 

151 

-225 

-59.8 

6/10 

243 

222 

-21 

-8.6 

6/23 

168 

113  | 

-55  | 

-32.7 

Average  Error  -22.6 


Average  Absolute  Error  36.0 
Standard  Error  36.9 
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TABLE  11 


Input  Data  for  Long-Record  Simulation 


Water  Years  1949-64 


Variable  | 

Interval 

Data  Source 

Precipitation 

Hourly 

Chicago  WB  City;  missing  data 
filled  in  from  Chicago  University. 

Potential  Eva- 
potranspirat ion 

Daily 

Double  duplicate  of  Chicago  O'Hare, 
1965-74. 

Temperature 

Daily 

Chicago  WB  City;  missing  data 
filled  in  from  Chicago  University. 

Wind  Speed 

Daily 

Double  duplicate  of  Chicago  O' Hare, 
1965-74. 

Radiation 

Daily 

Double  duplicate  of  Chicago  O'Hare, 
1965-74. 
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Simulated  Monthly  and  Annual  Flow  Volumes 
(Inches  of  runoff  for  28.5  sq.  mi.) 
Present  Land  Use  Conditions 
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TABLE  13 


Simulated  and  Recorded  Annual  Peak  Discharge 
West  Branch  DuPage  River  near  West  Chicago 


Simulated* 

Recorded 

WY 

Date 

Peak  Discharge 
(cfs) 

1  Date 

Peak  Discharge 
(cfs) 

49 

438 

1 

1 

50 

522 

51 

5-10 

364 

52 

1-14 

460 

53 

3-14 

196 

54 

3-25 

421 

55 

10-10 

947 

56 

10-6 

403 

57 

7-12 

1432 

i 

58 

7-2 

340 

59 

4-28 

1 

450 

60 

1-12 

473 

61 

9-14 

756 

9-26 

450 

62 

7-2 

378 

3-19 

361 

63 

l 

7-19 

351 

4-30 

217 

64 

7-7 

214 

7-19 

201 

65 

1 

7-14 

312 

2-7 

250 

66 

5-12 

1 

508 

5-12 

537 

67 

6-10 

1025 

6-10 

805 

68 

8-17 

739 

8-17 

340 

69 

6-9 

389  | 

6-9 

399 

70 

6-2 

438 

6-2 

521 

71 

12-11 

208 

2-20 

333 

72 

8-26 

687 

8-26 

715 

73 

12-30 

444 

12-30 

535 

74 

1-27 

314 

1-27 

522 

75 

4-19 

537 

76 

3-5 

557 

77 

3-29 

158 

78 

9-18 

468* 

♦  Constant  "present"  land  use  conditions  f  ^preliminary 
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Figure  t 

FLOWCHART  FOR  HVDROCOMP'S  HYDROLOGIC 
SIMULATION  PROGRAM  (HSP) 


SIMULATED  ACCUMULATIVE  RUNOFF  VOLUME  (inchat) 


Figure  3 

MASS  BALANCE  CURVE 

CALIBRATION  PERIOD 
OCT  64  -  SEP  69 
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Figure  7.2 
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Appendix  A 
Hypothesis  Testing 


Refer  to  Figure  7.1 

2  2 

(assume  =  o2  ) 


Observed 
X  =  2.6056 

Sx  =  0.1935 


Simulated 
X2  =  2.6553 

S9  =  0.2067 


n  =  18 


t  = 


(Xl-X2) 


n 


2 


26 


where 


'>  o 

sw 


and 


S“  =  { (n  -1 )  S  +  (n2-U  S2  ]  /  ( n 


S2  =  [ ( 18-1 ) (0. 1935)2  +  ( 26—1 ) (0.2067 )2 ]/( 18+26-2)  =  0.0406 


Sxrx2  = 


(0.0406  +  (0.0106),  55  =  0>0618 


18 


26 


,  (2.6056  -  2.6553) 
(0.0618) 


lor  a  -  0.05,  1  (  i-a/2)  (  n  t  hi  ,,-2)  ’(.976)(42) 


+n2-2) 


-2.018  <  -0.8046  <  2.018  therefore,  do  not  reject  II 


Refer  to  Figure  7.1  (Cont.) 


for  a  -  0.40,  >80) (42) 


0.851  Still  cannot  reject  H 


would  require  a  =  0.50,  t(.75)(42)  =  °-681  t0  re-ject  11 

2  2 
Hi  o  i  o 2 


F 


(0.1935)" 
(0.2067 ) 


0.0374 

0.0427 


0.8759 


for  a  =  0.05,  F(a/o)(n  lf  »  -1 )  =  F( .025) (17,25)  °‘3'° 


'(l-a/2)(n1-l,n,,-l)  "  ( .  975)  ( 17 ,25) 


:.44 


Note:  near¬ 

est  df  in 
table  are 
(15,24) 


0.370  <  0.8759  <  2.44  therefore,  do  not  reject  H 


for  a  --  0.50,  F(  .  25 )  (  1 7 , 25 ) 


0.712 


F(  ,75) ( 17 ,25) 


1.35 


0.712  <  0.8759  <  1.35  still  cannot  reject  H 
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Refer  to  Figure  7.2 


H: 


Mi  -  M2  (assume  a,  +  a2  ) 


Observed 


Xx  =  2.6706 


Simulated 


X2  =  2.6276 


Sx  =  0.1566 


S2  =  0.2101 


=  10 


n2  =  10 


t  '  = 


X1  -  X2 


2  2  V 

S1  S9  2 

(—  +  JL  ) 

ni  n2 


reject  if  t'  < 


“  (Vl  +  w2t2)/(wi  +  w2> 

or  t’  >  (w1t1  +  w2t2)/(Wl  +  w2) 


where  w  =  -J* 
1  ni 


tl  t( l-a/2 ) ( n  t  —  1 ) 


2  n. 


2  *( l-a/2 ) ( n  -1 ) 


Sample  X  S2 

1  2.6706  0.0245 

2  2.6276  0.0441 


n 

JO 

10 


w  = _ 

n 

0.0025 

0.0044 


for  a  0.05 
t 

2.262 

2.262 
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Refer  to  Figure  7.2  (cont) 


(2.6706  -  2.6276) 

t'  =  2  y  =  0.5189 

(0.1566)  +  (0.2101)  2 

[ -  ] 

10 


Vl  *  W2t2)  (0.0025X2.262)  +  (  0. 0044  )  (  2 . 262  ) 

(wi  +  w2)  “  (0.0025  +  0.0044) 


-  2.262  <  0.5189  <  2.262  therefore,  do  not  reject  H 


for  ol  =  0.50, 


t( l-a/2) ( 9) 


0.703 


-  0.703  <  0.5189  <  0.703  still  cannot  reject  H 


2  2 

Test  same  hypothesis  (H:  =  p9)  ,  now  assuming 

t  _  <x!  -  x2)  where  sf.  =  - —  +  - — 

-  s-  -  Xr5r2  ni  n2 

X1  X2 

and  S2  =  [  ( n^-1 )  S^2  +•  (n2~l)  S22]  /  (n^  +  n2  -  2) 


but  since  n^  =  n<;>,  the  above  simplify  to 


Si  -  ,  Ifillfsl’ 

xrx2 


s-  - 

xrx2 


[(0.1 566 ) 2  +  ( 0 . 2 1 0 1 ) 2 ] 
10 


0.0069 


Refer  to  Figure  7.2  (cont.) 


(2.6706  -  2.6276) 

y 

(0.0069)  2 


0.5189 


for  a  =  0.05, 


t( 1-a/ 2 ) {n^+n^-Z ) 


2.101 


-2.101  <  0.5189  <  2.101  therefore,  do  not  reject  H 
even  with  a  =  0.50,  still  cannot  reject  H 


Refer  to  Figure  7.2 


„  2  2 
H:  °i  =  °2 


F  = 


*1  (0.1566)  0.0245 


S0  (0.2101) 


2  0.0441 


=  0.5556 


for  a  =  0.05,  F,  /oW  =  Fr  no_Wo  =  0.248 

(a/2)(n1-l,  n  -1)  (.025)(9,9) 


F(l-a/2)(n1-l,  n2-l)  "  F(.975)(9,9)  "  4,03 


0.248  <  0.5556  <  4.03  therefore,  do  not  reject  H 


for  a  =  0.20,  0.410  <  0.5556  <  2.44  still  cannot  reject  H 
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Refer  to  Figure  7.3 


U1  ^2 


2  2 
(assume  =  o ^  ) 


Simulated  (26) 
X  =  2.6553 

S1  =  0.2067 


Simulated  (14) 
X2  =  2.6387 

S2  =  0.2050 


=26  n2  =  14 

(Xi  -  X  ) 

t  =  — -  where  S—  —  defined  as  before 


S2  =  [ (25) (0.2067)2  +  ( 13) (0.2050)2]  /  (38) 


—X 

X1  x2 


0.0425  0.0425  %  /z 

( — 26 —  +  ~T7—  >  =  °*°683 


14 


_  (2.6553  -  2.6387) 

t  “  0.0683  *  ^ 


0.0425 


for  a  =  0.05,  t,.  /oW  oX  =  t,  =  2.025 

( l-a/2) (n1+n2“2'  ( . 975 ) (38) 


-2.025  <  0.2429  <  2.025  therefore,  do  not  reject  H 


for  a  -  0.50,  t(>75)(38)  -  0.681 


-0.681  <  0.2429  <  0.681  still  cannot  reject  H 
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Refer  to  Figure  7.3 


H: 


2 

°2 


c  2 

^1 _  (0.2067) 

S22  (0.2050)2 


1.0167 


for  a  =  0.05, 

F(a/2)(n  -1,  n2-l)  =  F( .025 ) (25,13) 

F(l-a/2)(n  -1,  n2-l)  =  F( . 975) (25, 13) 

0.331  <  1.0167  <  2.54  therefore,  do  not  reject  H 


for  a  =  0.50 


F 


F 


0 


( .25) (25,13) 


*  0.684 


(  .75) (25,13) 


=  1.36 


.684  <  1.0167  <  1.36 


still  cannot  reject  H 


Note:  nearest 

df  in  tables 
are  (24,  12) 
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